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We analyze the structure of the high multiphcity events observed by the CMS collaboration at the 
LHC. We argue that the bulk of the observed correlations is due to the production of a pair of jets 
with pt > 15 GeV/c. We also suggest that high multiplicity events are due to a combination of three 
effects: high underlying multiplicity for collisions at small impact parameters, upward fluctuations of 
the gluon density in the colliding protons, and production of hadrons in the fragmentation of dijets. 
The data analysis is suggested which may clarify the underlying dynamics of the high multiplicity 
events and probe fluctuations of the gluon field as a function of x. 
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I. INTRODUCTION 

Recently the CMS collaboration reported the first data on the pp collisions at 7 TeV taken with a 
special high multiplicity (HM) triggerjll. The trigger selected events with multiplicity N^h > HO, in the 
pseudorapidity interval jryj < 2.4 and pt > 0.4GeV/c. This corresponds to the multiplicities ^ 7 times 
larger than in the minimal bias inelastic events {NminMas ^ 15). The reported studied the two-particle 
correlation in pseudorapidity, Arj — Tji — •q2, and in the azimuthal angle At/) = ipi — 4>2- The main focus 
of the discussion was on the observation of the same side positive correlation (ridge) at A0 ^ in a 
wide rapidity interval of 2 < |Ary| < 4. A similar correlation was observed previously in the heavy ion 
collisions at RHIC. 

The CMS ridge effect was discussed in a number of theoretical papers. However little attention was 
payed to a number of other remarkable features of the data. The aim of this note is to fill this gap and 
to suggest possible ways of analyzing the data which may clarify the dynamics of the HM events and 
provide a tool for studying fluctuations of the nucleon gluon density. 

Let us first summarize the basic observations of the CMS study. 

• The trigger selects events with average multiplicity which is approximately seven times higher than 
for the minimum bias inelastic pp collisions. Probability of the selected events is very small, 

Phm ~ 10-^ - 10-^ (1) 

• There is a very strong and localized in A77 positive two-particle correlation for /S.(f> ^ if two 
particles with 1 < < 3 GeV/c are selected. 

• For the same pt cut there is a strong positive two-particle correlation for At/) ~ tt for a broad range 
of A77. 

• Total correlated hadron multiplicity of hadrons in the same side 2 < jAT^j < 4.8 ridge is ^ 0.04 
(0.02) if both selected hadrons have comparable momenta in the range 1 < pt < 2(2 < pt < Z) 
GeV/c (the off-diagonal correlations were not reported yet). 

II. ORIGIN OF THE AWAY SIDE RIDGE AND CORRELATIONS IN A?7 ~ 0, A</) ~ 

The inspection of the correlation plot (Fig. 7 of ref. [ll) indicates that the total excess transverse 
momentum in the A0 ^ vr region is 



(2) 
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It is not possible to estimate in the same way the same side correlated transverse momentum since the 
correlation function for A0 = is off the vertical scale of the plot. However it is clear that the total 
excess transverse momentum in the A77 ^ 0, Aip ^ region is comparable to pb"-'-°-nce (ji]q^ 

Thus it appears that the bulk of the observed A77, Acf) correlations is due to production of two back 
to back jets with pt >15 GeV/c. This would explain both the strong narrow correlation in the towards 
region and a broad correlation in the away region. The later is due to a broad distribution over xi, X2 in 
the hard parton collisions at typical Xi ~ 2pt/y/s ^ 5 • 10^^. A priori this does not exclude a possibility 
that these jets are softer than the in-vacuum gluon jets observed say in the e~^e~' annihilation due to 
much higher gluon densities characteristic for these collisions (see discussion below) . This may be similar 
to the softening of the jets observed in the heavy ion collisions. Obviously, a softening of the jets would 
make it rather difficult to extract jets with such moderate transverse momenta from the data using 
existing jet finding algorithms. 

A natural question here is why pb"-'-°-nce jg large. Naively the production of jets with smaller 
Pt produced in the perturbative QCD (pQCD) regime with a higher rate should dominate since the 
centrality for jets with different pt is practically the same while the hadron multiplicity in the dijet 
fragmentation is a rather weak function of pt . It is possible that a large value of is related to 

the observation of the analysis ^ that pt at which pQCD starts to dominate grows with ^/s and may 
be as large as 8 GeV/c in the generic y/s— 7 GeV collisions. Though no specific dynamic explanation 
of this pattern exists so far, it appears likely that the pattern is related to the increase of the gluon 
density in the pp collisions with increase of y/s. If so, the effect should be stronger for collisions at small 
impact parameters which as we argue below dominate the HM collisions since in this case the average 
gluon densities of the partons involved in the collisions are significantly higher than in generic inelastic 
collisions. 

For the typical A77 ^ 2 we can estimate the invariant mass of the produced dijet system 

Mjet^ jet, = 2pt exp(A7//2) - lOOGeV. (3) 

The multiplicities of hadrons in the fragmentation of the gluon jets were studied in a number of 
experimental papers, see jsj and references therein - and found to be in good agreement with the pQCD 
expectations [J]. Using these data we can estimate that the charge hadron multiplicity in these events 
due to the gluon dijet fragmentation is ^ 30, and hence gives a substantial 25%) contribution to the 
total multiplicity. 

Presence of jets with large pt in the CMS HM events may be of relevance for the interpretation of the 
same side large | A77I ridge. Indeed, pQCD predicts a relatively small suppression of the particle density 



(integrated over particle energies) emitted in-between jets. The suppression is proportional to y^VcO^/Tr. 
Also, pQCD leads to the difference between in- and out-of plane emissions (a sort of "string effect"), in 
favor of in-plane emission - a ridge like structure. The actual numerical value is kinematics and color 



4 

flow dependent (for extended discussion of tlie QCD colierence plienomena in the jet production see Q]). 

It would be straightforward to find out how important the production of the back to back jets is in the 
HM events with the data at hand. This would allow also to check how much the pQCD "string effect" 
contributes to the same side ridge (which constitutes ^ 0.04 charged particle per event ioi 1 < pt < 2 
GeV/c). In particular one would be able to test predicted by pQCD dependence of the ridge on the 
the value of rji — ri2- Comparison with the ridge structure in the minimal bias dijet production for 
Pt > 20GeV/c would be useful as well. 



III. DYNAMIC MECHANISMS FOR GENERATING HIGH MULTIPLICITY AND HIGH pt 

JETS 



It is natural to expect that the average hadron multiplicity should monotonously increase with a 
decrease of the impact parameter, b. Hence the HM trigger of the CMS should correspond to the 
collisions at very small b. We can use the analysis of Q to estimate the probability that an inelastic 
collision occurs at small impact parameters. We find that the probability for b < 0.2 fm is small ~ 2% 
but still much larger than Phm (EqllJ. 

Let us now estimate the average multiplicity for collisions at 6 ^ for ^/s = 7 TeV. The inclusive dijet 
trigger selects collisions at a median b « 0.6 fm The underlying multiplicity for such collisions is 
about a factor of two larger than in the minimal bias non-diffractive events, see discussion and references 
in 0,0]. So for ^ the enhancement should be somewhat larger. 

A rough estimate can be made using information about 5-distribution of minimal bias events for which 
median & is ^ 2 than for the dijet events [sl. Experiments excluded diffractive events when calculating 
minimal bias hadron multiplicity. The diffractive processes which constitute about 30% of the total 
inelastic cross section mostly occur at large impact parameters. Taking this into account and comparing 
the b distribution in non-diffractive minimal bias events and in dijet events one can " subtract" the large 
b tail of the dijet distribution. As a result we estimate that the average multiplicity for collisions with 
b < 0.6/to is about 2.5 times higher than in minimal bias non-diffractive events. One may expect a 
further increase for 6^0. 

Hence we conclude that a scenario where small Phm is solely due to the small probability of small b 
collisions and a small probability of selecting dijets with high pt would lead to a large multiplicity on 
the scale of 

Nch — Njetjet + Nunderlying > 70. (4) 

In this estimate we have assumed for illustration that dijet multiplicity is of the order one. However the 
obtained Nch is still significantly smaller than the observed one. Thus HM trigger selects the tail of the 
distribution of the multiplicity for central collisions. 
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It is natural to expect that to generate a higher muhiphcity one needs to take into account fluctuations 
of the strength of the gluon field on an event by event basis. It was demonstrated in [7] that one can 
relate the dispersion of the strength of the gluon field at small x to the ratio of inelastic and elastic 
vector meson production at t = 0: 



Wn = 



(G)2 



dt / dt 



(5) 

t=o 



A model of global fluctuations was proposed in [7| which allowed the explanation of the magnitude 
of the experimental ratio in Eq. [S] It took into account the QCD DGLAP evolution of ujg. For the 
discussed CMS kinematics the model jj| leads to ujg ^ 0.1. 

We can define the probability for a fluctuation in a nucleon to have the strength of the gluon held > r 
times larger than the gluon parton density which is the average over all configurations in the nucleon as: 

P{r)^ J d<j^^^^0iGix,Q'\a)-rG{x,Q')), (6) 

where a labels different configurations in nucleons. 

Making a simplifying assumption that the distribution over the strength of the gluon field is Gaussian 
we can estimate 

P(1.5) - 3%. (7) 

However it is likely that the probability of the large fiuctuations is larger than given by Eq. [T] Indeed 
the analysis ^ of the data on production of four jets (7 + 3 jets) using information about general- 
ized parton distributions indicates presence of the positive multiparton correlations in nucleons. Such 
correlations are likely to increase P{r > 1.5). 

Hence an overall probability for the interaction of two protons to occur at very small impact parameters 
b < 0.2 fm and with both colliding nucleons in configurations having gluon density > 1.5 times larger 
than average, should be > 10~^ which is comparable to the frequency of the events selected by the HM 
trigger (Eq. [1]). This scenario would yield a larger value of Nch than the estimate of Eq. |4] due to a 
higher gluon density in the colliding configurations both due to soft and semi-hard interactions. Also it 
would increase the rate of binary parton collisions per event. 

To estimate the rate of dijet production in 2— > 2 collisions we consider first the case of inclusive dijet 
production. The normalized distribution over b of the inclusive dijet production, P2{b) was calculated 
m through the convolution of the generalized gluon parton distributions which are measured in 

the exclusive hard processes at HERA. For the dipole parameterization of the two gluon form factor 
iF,git) = il-t/mlix))~'): 

P.ib) ^ ^ Ksirugb). (8) 
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For the central collisions where 6 ^ it is straightforward to find the ratio of the multiplicities in the 
6 ~ collisions and in the minimal bias collisions: 

Ro = P2{0)<jUpp) - -^<7„App). (9) 

iZTT 

For a; ~ 5 • 10~^, Mg w IGeV^, leading to a significant enhancement of the dijet rate R — 4.6. 

There are two types of fluctuations which modify R - one is fluctuations of g^ix, Q'^\a), and another 
is fluctuations of the size of the overlap area, S, of the colliding protons. For example, let us consider 
a model of a nucleon where the gluon field is mostly in the plane of three valence quarks - " a pancake 
shape" . In this model S for the collisions when the pancakes collide edge on edge is much smaller than 
in average. If we consider collisions in which the three quarks in each of the nucleons are aligned along 
the reaction axis S would be even smaller. Since the rate of the jet production is cx l/S* such collisions 
would lead to a large enhancement of the jet multiplicity. 

Inclusion of the fluctuations would increase i? by a factor ri ■ r2'. 



R = i?orir2, (10) 

suggesting that the enhancement could be as large as 10 for HM events, where 

^^^^ ^ gN{xi,Q'^\a)giN{x2,Q'^\cr) (S) ^^^^ 

gN{xi,Q'^)giN{x2,Q'^) S 
The multiparton 4— > 4 hard collisions occur at even smaller impact parameters than the 2— > 2 hard 
collisions. Consequently, the enhancement of the four jet production due to the 4— > 4 interactions in 
HM events as compared to the minimal bias events should be larger than in the 2-^2 case. We can 
estimate the enhancement factor in the independent particle approximation using expression for Pi{b) 
from Including the effect of fluctuations in the gluon density we find 

7 

Rijet = 77^™o^m''l''2f3?'4, (12) 
OOTT " 

which is several times larger than for the two jet case. However the independent parton approxima- 
tion underestimates the rate of four jet production [2], so the enhancement factor may depend on the 
mechanism of parton-parton correlations and differ ffrom the square of the enhancement for the dijet 
production. 



IV. TESTING COLOR FLUCTUATION CONJECTURE 

Above we have made a conjecture that the HM trigger selects collisions of nucleons in configurations 
with larger than average gluon density. In other words we suggest that HM and high gluon density 
fiuctuations go hand in hand. 
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It is rather difficult to clieck this conjecture using the data analysis procedure of ref. [1] since one 
could always question whether the HM trigger by itself was selecting events with extra hard collisions. 

Here we discuss several possible strategies which address a possible dependence of the gluon fluctu- 
ations on X, virtuality of the gluon and overall transverse size of the projectiles. Note here that so far 
the gluon fluctuations were modeled only in [7]. This model focused on the small x region and led to a 
decrease of co^ with increase of x. In this model significant cancelations are present for the enhancement 
factor given by Eq. [TT] since in this model the upward fluctuation of the gluon field are related to the up- 
ward fiuctuations of the overall size for configurations. However the model does not include fiuctuations 
of the shape which may not change the gluon densities, for example "pancake" or "cigar" fluctuations 
discussed in the previous section. Such effects may lead to significant fiuctuations of gig2/S for a wide 
range of x. 

One option to probe gluon field fluctuations is to look for dijets at the central rapidities in matching 
back to back narrow cones but to use in the HM trigger only hadrons outside these cones. This way one 
would be able to determine how much gluon flelds in two nucleons are enhanced in the overlapping area 
as compared to the scenario where no fluctuations are present. 

Another option is to consider the multiplicity of the dijet production, N2j in the region |?7| > 3.5 ^ 3.0 
as a function of the multiplicity at central rapidities, |r/| < 2.4 - Ncentrai- The gap between the two 
regions would minimize the leaking of hadrons from hard collision to the central region. This problem 
is of importance only if rji > 0, 772 < 0. Both the CMS and the ATLAS can measure jets in this 
kinematics using forward calorimeter components of their detectors. In this kinematics one would test 
the fluctuations of gluons either with say xi ~ 10^"^ and X2 ^ xi (same side jets) or if xi ^ X2 > 0.05 
(opposite side jets). 

One can expect that N2j (Ncentrai) would grow monotonously with increase of Ncentrai due to increase 
of the centrality of the collisions. There should be a substantial increase of N2j (Ncentrai) between 
Ncentrai = Nmin.bias and Ncentrai = '^N^inMas Corresponding to typical underlying multiplicity for dijet 
events (which reflects a more central nature of the hard collisions) . If the gluon fluctuations would not 
kick in, N2j would saturate at the value 

N2j = XN2j(N„,,nMas), (13) 

with A ^ 4, cf. EqE] with a small residual Xi dependence due to the broadening of the transverse 
distribution of gluons with decrease of Xi 0, Q 

Fluctuations may lead to a continued increase of N2j for Ncentrai 3> '2Nmin.bias due to the fluctuations 
of the gluon field. The rate of the increase may depend on Xi and pt of the jets. In particular, in the 
model of Q the increase of N2j for b ^ events should be weaker for 771 ^ —772 ^ 3 than for r/i ^ 772 ^ 0. 

In principle, by combining measurements in three discussed kinematic ranges, it may be possible to 
disentangle the x-dependence of the gluon fiuctuations. 
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V. CONCLUSIONS 

In the presented scenario HM events occur when the protons interact at very small impact parameter 
in special configurations with enhanced gluon field density. Production of two back to back jets with 
Pt > 15 GeV/c in these events is strongly enhanced as compared to the minimal bias events. 

To test this hypothesis it is necessary to study the rate of jet production in the HM events in the 
vicinity of the test particle, and the degree of correlation of such jets with a recoil jet production. 
Measurement of the ratio of the probability of high pt jets in the HM at a wide range of rapidities and 
minimal bias events would be especially revealing and may test the conjecture of the gluon fluctuations 
in nucleons and study x-dependence of the strength of these fluctuations. 

To check the relevance of the pQCD string effect for the explanation of the same side ridge it is 
necessary to perform a comparison of the HM correlations with the correlations in pp scattering with 
production of dijets with pt > 20 GeV/c. Alternatively the ridge could be due to a selection of rare 
configurations in the protons of a complicated transverse shape. 
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Note edited in proof. Very recently the ALICE collaboration reported preliminary data on the depen- 
dence of the J/ip multiplcity in pp collisions at ^/s — 7TeV on the hadron multiplicity at rj close to zero. 
For the highest measured hadron multiplicity which is about a factor of 5 larger than the minimal bias 
multiplicity they report an increase of the J ftp multiplicity of about 4 5. This is close to the estimate 
of Eq.(9). 
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